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Abstract 

Jet cross sections in deeply inelastic scattering in the case of transverse photon exchange 
for the production of (1+1) and (2+1) jets are calculated in next-to-leading order QCD (here 
the '+1' stands for the target remnant jet, which is included in the jet definition for reasons 
that will become clear in the main text). The jet definition scheme is based on a modified 
JADE cluster algorithm. The calculation of the (2+1) jet cross section is described in detail. 
Results for the virtual corrections as well as for the real initial- and final state corrections 
are given explicitly. Numerical results are stated for jet cross sections as well as for the ratio 
cr (2+l) jet/ CT tot th at can be expected at E665 and HERA. Furthermore the scale ambiguity 
of the calculated jet cross sections is studied and different parton density parametrizations 
are compared. 
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1 Introduction 



Recent results from E665 and HERA show that events with a clear jet-like structure are present in 
deeply inelastic electron proton scattering j|, With sufficient luminosity it should therefore 
be possible to study jet cross sections to use them for a test of QCD and an independent 
determination of its fundamental parameter Aq^rj- Jet cross sections may even be useful to 
extract some information on the gluon density at small x, because the gluon density is important 
in (2+l)EJ jet production. 

Because of the strong scale dependence of fixed order cross sections calculated in perturbative 
QCD, the possibility of a large size of the corrections and because a determination of a s must be 
based on a next-to-leading (NLO) order cross section, the calculation of higher order corrections 
is well motivated. Since now experimental results are available, it is worthwhile to give a detailed 
account of the technical problems of the calculation of the (2+1) jet cross section. A second goal 
of the present work is to study the jet cross sections in detail numerically. In deeply inelastic 
scattering, the 0(a s ) corrections to the 0{cP s ) Born term are well known (see ||, |3|, ||, || [7J8], J|). 
In addition, the O(o^) Born terms for the production of (3+1) jets have been calculated [1C|, [ll]]. 
In this paper the calculation of the cross section for the production of (2+1) jets to 0{a 2 s ) in 
the case of transverse photon exchangee is described. For more details of the calculation see 
pLpfl. The contribut ions with a transverse exchanged photon dominate the cross section, as is 
shown in Section ||[ More recently, the other parity-conserving helicity cross sections have been 
calculated based on p2| , see [14]. 



One of the main features of hadronic events in high-energy collisions is the pronounced jet- 
like structure. These jets are attributed to the production of partons in the fundamental QCD 
process p|. Due to the presence of infrared divergences in QCD, a suitable prescription has to 



be given in order to define the finite parts that arise from divergent terms after the singularities 
between real and virtual corrections have been cancelled. Such a prescription is related to the 
jet definition that is used on the parton level to calculate jet cross sections. In Section ^ jet 
definitions in deeply inelastic scattering are discussed and the jet definition scheme based on a 
modification of the JADE cluster algorithm that is used in this calculation is defined. 

In Section |H the calculation of the (1+1) jet cross section in next-to-leading order is reviewed 
and the results for the (2+1) jet Born terms are stated. In Section I] the results for the virtual 
corrections to (2+1) jet production are given. In Sections ||| and | the calculation of the real 
corrections to (2+1) jet production is described (separated according to singularities in the final 
and initial state). The sum of virtual and real corrections gives the finite jet cross section after 
renormalisation of the parton densities. The flavour factors are listed in Section [7[ In Section 
H the numerical results for jet cross sections and for the ratio °"(2+l) jet/ a tot are presented. 
The dependence of the jet cross sections on the renormalisation and factorization scales, the 
dependence on the jet definition scheme and the dependence on the parametrization of the 
parton densities is also studied. The appendix contains explicit results for massless 1-loop 
tensor structure integrals, phase space integrals and the virtual and real corrections. 



1 The target remnant is counted as a jet, so "2+1" stands for the production of 2 partons in the hard QCD 
process, possibly accompanied by additional soft or collinear partons. 

2 Here, the cross section for transverse photon exchange is defined by the helicity cross section obtained by a 
contraction of the hadron tensor with the metric tensor (—g^)- 
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2 Jet Cross Sections 



Jets must be denned in terms of experimentally observable quantities. An experimental event 
is characterised by the energies and momenta of the outgoing particles. Since at high energies 
the multiplicity of the events is large and since presently there is no practical way (based on 
QCD) to describe hadron dynamics on the level of observable particles, one has to try to extract 
information from experimental data in a form that can be compared with theoretical results 
from perturbative QCD. At e + e~-colliders it has been observed that outgoing hadrons very 
often appear as clusters of particles concentrated in a small cone in momentum space. These 
clusters were called jets. Given an experimental event and a resolution parameter c (the "jet 
cut"), a suitable algorithm is applied to the event giving the number of jets of the event and 
the particles associated with each of the jets. Therefore, the algorithm that is used defines what 
is meant by a "jet". To compare experiment and theory, one must use the same algorithm in 
theoretical calculations. One should expect that experiment and theory are comparable as long 
as the same algorithm is used in the experimental analysis and in the theoretical calculation. 
Of course, the problem is that in realistic events the final state consists of hadrons, whereas 
in theoretical calculations (based on QCD) the outgoing "particles" are partons. Therefore the 
crucial hypothesis is that jets on the hadron level and jets on the parton level can be identified. 



The first jet definition that has been given is that of Sterman and Weinberg [15 1. It is 
based on cones in momentum space defined by an opening angle 5 and an energy fraction e. 
This definition, which is well suited for e + e~-annihilation for a small number of jets becomes 
complicated if a larger number of jets is produced. In addition, it is not Lorentz invariant (this 
is not a problem in the case of e + e~-annihilation, since here the CM system is a unique frame 
of reference). Later, another type of algorithm was proposed, the cluster algorithm first used 
by the JADE collaboration [16]. This algorithm combines succesively two particles into a jet, if 
their invariant mass squared is smaller than a fraction c of a typical mass scale : 



s i.j 



< cM 2 . (1) 



In e + e -annihilation, c is of the order of 10 2 , and M 2 is set to Q 2 , the total invariant mass of 
the hadronic final state. 

For hadron colliders, jet definitions in terms of "cones" in rapidity and azimuthal angle are 
favoured (UAl-type algorithms). Such a jet definition singles out a particular axis (namely, the 
direction of the two colliding beams). In the case of pp-events the incoming partons are assumed 
to have only small transverse momentum, and therefore the situation is symmetric with respect 
to this particular axis. 

The situation is quite different for electron-proton colliders. Here the interaction is mediated 
by the exchange of a virtual photon with momentum q (and Q 2 := — q 2 > 0). This photon hits 
a proton with momentum P. Therefore, the interaction should be described in the CM system 
with P + q = 0. This system varies from event to event, and this is the reason why a Lorentz 
invariant jet definition should be used [|ll]] . 

Before a suggestion for a suitable jet algorithm in eP-scattering is given, one should have a 
closer look at the target remnant. An interesting question is: "Is it reasonable to include the 
target jet in a jet analysis?". One should consider the process in fig. [l]. This Feynman graph 
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describes initial state radiation of a gluon with momentum p2 from the initial quark line with 
momentum pq. It is assumed that the gluon is emitted collinearly with a large energy in the 
direction of the incoming quark causing a strong enhancement of the cross section for this process 
because of the pole of the quark propagator. Since all partons are assumed to be massless the 
gluon will go in the same direction as the target remnant. Under the assumption that "parton 
jets" roughly correspond to "hadron jets", there will be no possibility to disentangle the hadrons 
from the debris of the proton and those coming from the fragmentation of the gluon. To be 
consistent, one therefore should define "parton jets" in the following way: If an outgoing gluon 
can be separated from the remnant by a suitable condition (e.g. invariant mass), the gluon and 
the remnant are considered to be two jets with momenta P2 and p r . If the gluon and the remnant 
cannot be separated, they count as one jet whose momentum p* = P2 + Pr is the sum of the 
gluon momentum and the momentum of the remnant. In an experiment, however, one cannot 
measure the momentum of the remnant directly since most of the hadrons from this jet are lost 
in the beam pipe. The momentum of the target remnant jet must therefore be determined in 
an indirect way. The jet algorithm should also have the property that all collinear singularities 
are treated in such a way that they factorize the corresponding Born term. This allows for a 
process independent definition of the renormalised scale dependent parton densities. 

For the experimental analysis it is proposed to use a modified JADE cluster algorithm 
(mJADE algorithm) consisting of two steps (see also JTl|): 

(1) Define a precluster of longitudinal momentum (in the direction of the beam pipe) p r that is 
given by the missing longitudinal momentum of the event. 

(2) Apply the JADE cluster algorithm to the set of momenta 

{pi,P2,---,Pn,Pr}, (2) 

where pi,P2, ■ ■ ■ ,Pn are the momenta of the visible hadrons in the detector and p r is the mo- 
mentum of the precluster. 

It remains to define the order of magnitude of the jet cut and the mass scale M 2 to be used 
in the mJADE algorithm. The jet cut should be such that it is small enough to ensure properly 
separated jets in the detector (otherwise the jets coming from the real corrections could become 
too broad), but large enough to avoid large logarithms that could spoil a fixed order result in 
perturbation theory. One can think of cM 2 as a new mass scale that is needed to specify the 
cross section completely. In a deeply inelastic process there are several mass scales given by 
an event, namely the virtuality Q 2 of the exchanged vector boson, the invariant mass W 2 of 
the hadronic final state and, in general, several invariant masses p 2 ± from transverse momenta 
of outgoing particles. Therefore large logarithms are expected if the quotient of any two of 
these scales becomes small (or large). So one should avoid kinematical regions where this could 
happen. Because of the jet cut an additional scale cM 2 enters the calculation, where one could, 
for example, use Q 2 , W 2 or some p 2 ± as M 2 . Since it is suggested to include the target remnant 
in the jet analysis, it is natural to use the scale M 2 = W 2 in the jet algorithm (analogous to 
the situation in e + e~-annihilation). If the scales pf ± are omitted from the discussion, the scales 
Q 2 = S H yx B , W 2 = S H y(l — x B ) and cW 2 = cS H y(l — x B ) are relevant. Here \fS^ is the total CM 
energy of the collider, x B = Q 2 /2Pq is the Bjorken variable and y is the usual lepton variable 
y = Pq/Pk, with P the proton momentum and k the momentum of the incoming electron. 

This is not the only reasonable choice for the mass scale used in the the jet definition. In 

fact, in Section [3| it is pointed out that a scale like Q 2 or (w a Q^^S H y 1 01 /3 )) with parameters 
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a, /3 may be reasonable if the parton densities have to be probed at small x. In contrast, the 
considerations here are based on the simple observation that the observable invariants (including 
those with the remnant jet) sum up to W 2 , and it is somehow natural to use this particular 
scale. For another possible jet definition in eP scattering see [jlTj. 

In the cross section large logarithms of Q 2 /W 2 « x B , cW 2 /W 2 = c and cW 2 /Q 2 « c/x B (the 
approximation of small x B is made since most events are expected in this region) are expected. 
For cuts of the order of 0.01 the logarithms in c should be comparable to those encountered in 
e + e~-annihilation because of a similar structure of the matrix elements. 

Berger and Nisius have studied the effect of the inclusion of the target remnant in the jet 



analysis [18] by using the Monte Carlo generator LEPT05.2. They come to the conclusion 
that the correlation of the number of jets on the parton level and the number of jets after 
fragmentation is much stronger if the remnant jet is included. For more details, see [11]. 



3 (1+1) and (2+1) Jets: Cross Sections to G(a s ) 

In this section the calculation of the (1+1) jet cross section in NLO is reviewed. This is done 
for two reasons: these results are needed to calculate jet rates, and their calculation serves as an 
illustration of the more complicated case of the NLO corrections to (2+1) jet production. As a 
byproduct, one gets the results for the (2+1) jet Born terms. The total cross section to 0(a s ) 
for photon exchange has been calculated in and the (1+1) jet cross section to this order for 
all neutral and charged current processes can be found in ||. Here the special case of the (1+1) 
jet cross section for the exchange of a transverse photon is discussed in detail and results for a 
longitudinally polarised photon are stated. It is safe to focus on the QCD corrections for the 
transverse polarization of the virtual photon, since the cross section for longitudinally polarized 
photons contributes only about 20% to the Born term cross section (see Section ||), and this is 
expected to be true for the relative contribution of the longitudinal terms to the NLO corrections 
as well. 

The cross section for eP-scattering differential in x B and y is given by 

da H _ ^^dt f . c UKy(S H x B )-t{y{l-y))^^ 2 l 1 



dx B dy ^L r . Z J ' iU r(l-e) ~ 2S H x B 

■ ^ S parton ^Jfe^M Q2 JfF j dn ' ^ H ^ ( 3 ) 

where d = (4 — 2e) is the space-time dimension (e / regularises the ultraviolet and infrared 
divergences fl9|, ^0[ |2l|), pi is a mass scale for making the coupling constants dimensionless in d 
dimensions, rfPSp^ rton is the n-parton phase space 

barton = (2*)" J[ ^S** to, + *-£»), (4) 

i=l ' i=l 

Q' is the volume element of d — 3 angles specifying the direction of the outgoing lepton relative 
to the outgoing partons, and 

N' = f dn' (5) 



Q 
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is a normalisation constant. The incoming parton carries a fraction £ of the proton momentum, 
the are the bare parton densities for partons of flavour i, 

= k^k' u + k u k'» - ktfgT (6) 

is the lepton tensor and H^ v is the hadron tensor (including coupling constants, colour factors, 
etc.). 

The integration over fi' can be performed. It can be shown by a direct evaluation of the 
integrals in d-dimensional space that 

1 1 f^nu, _ i + (i-y?-ey 2 



Q 2 N'J 1 ~ 2(1 -e)y 2 [ 9 € " ' 



, 4(1 -e)(l-y) + ! + (!- y) 2 - ey 2 ( ^ 

+ 2(1 -e)y 2 {e R +£ h (7) 



where 



1 



q Q 
1 



Q 2 

Ax 



2 



£ - -Q2P0PO1 (°) 

and x p := x B /£. Because of current conservation q^H^ = one has 

e^H^ = e^H^ = 0. (9) 
With the definition trH := h g := (—g^ v ) H^, ho := e^H^ one obtains 

l + (l-y) 2 - g y 2 A{l-e){l-y) + l + (l-y) 2 -ey\ 
= 2(1 - e)y 2 h » + 2(1 - e)y 2 ^ {W) 

'i d£ r _ (AttY (S h x b )-^ (y(l - y))- ^ 2 1 1 



By defining 



Q: 



x s £ ' ' r(l-e) 2fe 

" dPS parton e 2^2e( 27r )2d ^ A t 11 ) 



with A G {y, 0} one arrives at 

da H 1 + (1 - y) 2 - ey 2 



cfe s dy 2(1 - e)y 

4(1 - e)(l - y) + 1 + ( 1 - y) 2 - ey 2 
2(1 - e)y 



t" ^ 1^1 a o- (12) 
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In the literature cross sections ajj, ai for unpolarized and longitudinally polarized photons have 
been defined. They are related to the definitions used in this paper by a g = 2(1 — e)ajj — <jl, do = 

Now the special case of (1+1) jet production is considered. Fig. |2| depicts the Feynman 
diagram to O(a s ). The dia gram for the 1-loop virtual correction is shown in fig. |3[ The result 
for the sum of these diagrams is well known : 



da, 



Born&virt. 



dx B dy 



dx 



E J ?r 2,M 

• J x B Xp \ Xp 
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(4tt) £ (S H x B y* (y(l - ^ 2 1 1 



rfi 



-a - 



2 SfjXg 



l + (l-yf 



2(1 - e)y 2 



+ 



47T / U 2 



2vr4(l-e) Q 2 

; r(i- 



Q 2 I T(l - 2e) 2vr 



Cf — n 



+ ^^ F (-8-2C(2)) 



rfi 



T(l - 2e) 2vr e 



" " s -P q ^(x p )\ + 0(e). 



8(1 - x p ) 
(13) 



The integration over £ is rewritten in terms of the variable x p = x B /£, and the symbol 8 q i restricts 
the summation to quark initiated terms. ((2) = 7r 2 /6, and Qi is the charge of the quark with 
flavour i normalised to e. is the renormalisation scale and Mf is the factorization scale. The 
pole terms in e proportional to 6(1— x p ) are infrared singularities that will cancel against infrared 
singularities in the real corrections, and the term proportional to the Altarelli-Parisi splitting 
function P q ^i(x p ) will cancel against collinear singularities in the real corrections. Note that the 
bare parton densities are already expressed in terms of the renormalised parton densities 
0, |22] in the MS-scheme 



f! eU (e,M 2 ) 



du 



5 



6ij6(l - u) + 



a s 
2vr 



n T l - 2e) 



The Altarelli-Parisi kernels are 

1 + u 2 





-«(«) 


= C F 




-?(«) 




p 9 + 




= 2iV c 






1 r 

= 2i" 



l + (l-u) 2 



u 

1 



(1 



u 2 + (l 



+ - + «(1 
u 



u 



+ 



Nf is the number of quark flavours. 



u 
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A c - -JV> 6(1 - u) 



(14) 



(15) 



Now the real NLO corrections to (1+1) jet production will be calculated. The Born terms of 
0(a s ) have to be integrated over the phase space region that, by the jet definition, is considered 
to be a (1+1) jet region. In a first step, suitable variables are defined, then the (1+1) jet region 
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is determined and finally the integration is performed. The phase space for the production of 
2 partons is constructed in the following way. As usual, the variable 

z = iv (16) 

is defined, where pi is one of the outgoing partons. By defining t = su/W 2 (p2 is the momentum 
of the second outgoing parton) and a = x B + (1 — x B )t one obtains 

J (iPS (2) = J dPS* (2) ^(£-a), 



/ 



where the fact that £ = a if pq = £P is being used. By means of the factor — a) the integral 
/ /(£) m the cross section formula (||) can be performed trivially. The ranges of integration 
are z £ [0, 1] and t £ [0, 1]. The invariants Sjj = 2pipj expressed in terms of z and t are 

s i = S H y(x B + (1 - x B )t)z, 

S02 = S H y(x B + (1 - x B )t)(l - z), 

812 = S H y(l - x B )t. (18) 

The momentum of the target remnant is p r = (1 — £)P. The observable momenta are p±, p2 and 
p r . Invariants for these momenta are defined by Uij = 2pipj /W 2 . In terms of the variables z 
and t they read 

u r \ = (1 — t)z, 

U r2 = (l-t)(l-z), 

U12 = t. (19) 

Let us define the (2+1) jet region by the condition Sij > cW 2 , with i,j £ {1, 2, r}. It is easy 
to see that in the (t, z)-plane where the allowed phase space is the square [0, 1] x [0, 1] the (2+1) 
jet region is a triangle given by the following conditions: 

(a) c < t < 1 - 2c, 

(b) z c (t) < z < 1 — z c (t), where z c {t) = 

The region within the square surrounding this triangle is therefore the (1+1) jet region. All 
regions where the cross section becomes singular are within the (1+1) jet region thus allowing 
the factorization of the collinear singularities and their absorption into the renormalised parton 
densities. It should be noted that for (2+1) jet production the minimal momentum fraction £ of 
the parton densities that can be probed is £min( x s) = x B + (l — x B )c because of the cut condition 
on si2- The minimum is £min = c for x B — > and is therefore of order 0.01. If one wants to 
use the gluon initiated (2+1) jet events to determine gluon densities at small x a different scale 

for the jet definition could be chosen. If the scale is Q 2 or (w a Q^ '^S H y l a ^ , the minimal 
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£ depending on x B turns out to be £, m i n (x B ) = x B (l + c) and £ m i n (x s ) = x B + (1 — x B ) a x B @c, 
respectively, and for sufficiently small x B and c both give small £ m j n . If x s is fixed, these 
definitions are of course equivalent to the choice of a much smaller cut c in the jet definition 
scheme based on W 2 , and it is therefore important to check whether the (2+1) jets are still 
calculated in a regime where perturbation theory is still valid, i.e., whether the (2+1) jet rate 
is small enough compared to the total cross section. 

After this short digression the results for the Born cross sections for the production of 
2 partons are stated. The relevant diagrams are depicted in fig. || with the graph G replaced 
by the diagrams of fig. ||. Of course, in fig. [| one would have an additional diagram from an 
incoming antiquark, but since the resulting expressions are (except for the charges) the same as 
those for incoming quarks, these diagrams are omitted in the sequel. 

The traces of the hadron tensor for quark and gluon initiated processes are 

tr ^Born, q inc. = LlQ^i 1 ~ e )CpT q , 

tr -#Born, g inc. = ^lQj8(l - e)-T g (20) 



with 



'£02 £12 \ 2Q 2 Sqi 
.512 S02/ S02S12 



J-( ( l- e) f^ + f-)_^_ 2e \ ( 21 ) 



and 



L X = (2vr) 2 VeV 4e - (22) 



Po is always the momentum of the incoming parton, and p± is the momentum of the outgoing 
quark. 

These formulae are already averaged over the colour degree of freedom of the incoming 
partons. The factor of 1/(1 — e) in the expression for T g comes from the fact that a gluon 
has 2(1 — e) helicity states in d = (4 — 2e) dimensions compared to only 2 helicity states in 
4 dimensions. 

Now the integration of the cross sections over the (1+1) jet region of the 2-parton phase 
space can be done by subtracting the (2+1) jet cross section from the total cross section to 
0{a s ). The technical reason is that this subtraction makes the explicit universal factorization 
of the singular parts in the form of a product of an Altarelli-Parisi splitting function and the 
Born term more transparent than the direct integration over the (1+1) jet region. 

The total cross section from the real corrections for transverse photons and quark initiated 
processes is ]|]: 

A^° M ' q _ j2 f'^Efrenf^B ^ (^f {S H x B )^ {y(l - y))^ ^ 1 



Q: 



dx B dy . J x B Xp J l-Xl ^) 2 Sh x b 

1 + ( ]-^- ey2 2* "(i - e) Q\ (*£)' 2: c F s„ 

2(1 - e)y 2 y ' 1 \ Q 2 J T(l - 2e) 2tt q 
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-~5(1 -x p ) + - 



+2 



ln(l - x p ) 

1 Xn 



+3-x p + -<5(l 



-— Pq ^ q + 35(1- 
- (1 + x p ) ln(l 



+ 0{e). 



Xr, 



1 + x 



2(1 



— In x r 



(23) 



The "+" -prescriptions are defined in |23| (they are implicitly used for functions defined on the 
interval [0, 1]) and are the result of the subtraction in the collinear regime. The bare parton 
densities are expressed in terms of the renormalised ones and terms of 0(a 2 s ) are dropped. The 
corresponding expression for the gluon initiated process isi 



dx B dy 



2N f 

E 



1 dx 



-f 

J < 




(47r) e (S h x b )-* (y(l-y))-^ 



e " 4e 2 1 1 
a 



i + (i-y) 5 



2(1 -e)y 



2-R 



q^g 



- 2 2vr4(l-e)Q 2 ( 

(x p ) + ((1 - xpf + x: 



r(l - e) 2 S H x B 

'47T// 2 V r (!- e ) a s 1 1 



Q 2 



In- 



1 



T(l - 2e) 2?r 2 1 - e 



+ 0(e). 



(24) 



Note that the sum over quark flavours reflects the different flavours that are produced in this 
process s 

In the sum of (|l3|), ( p3[ ) and (p4|) the infrared and collinear singularities cancel, and one is 
left with a finite total cross section to 0{a s ). To obtain the (1+1) jet cross section one has to 
subtract the (2+1) jet cross section from the total cross section. Since there are no singularities 
in the (2+1) jet region, e can be set to 0. The integration of (|20| ) over the (2+1) jet region is 
not difficult. If all contributions are summed up one arrives at the finite (1+1) jet cross section. 
The quark initiated part is given by 

2N J rl An. . All l + (l-y) 2 



4 1+1),t ' q 



dx B dy 



E 



dx 



- q! fi 



rcn 



-,M 



/ 



a 



2 SiiXu 



2y? 



+2 



ln(l 



1 — x 



(-8-2C(2))<y(l-z p ) 
- (1 + Xp) ln(l 




•2vr-4 



1 + xi 



+3-x p + -5(1 -Xp) + — In j^P q ^ q (x p ) 



1 



1 1 



2 1 



+ 



x P + 2- 



PJ 2(1 

g_ 2 

/ 

1 - 2z c {t(x p ))) 
ln 1 - ^(sp))' 



1 



lnx r 



Xr, 



(25) 



Compared to |4| there is a difference in the finite parts. It results from the average over the helicities of an 
initial gluon. Here it is assumed that the gluon has d — 2 = 2(1 — e) polarization states instead of 2. 

4 To make the cancellation of the collinear divergence transparent, a "double counting" is included which is 
cancelled by a factor of 1/2. 
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and the gluon initiated processes contribute 



d4 1+1)>t ' g 

dx B dy 



2N f 



dx 



i—l Jx B x p y x p 



2 1 1 1 + (1 

Q -- 



2 S H x B 



2y 2 



1 «5 

2 2vr 



((l_ Xp )2 +a; 2) hn I 



-^-l") +2 1n^P f 



27T-4 



- (- (1 - 2 c (i(x p ))) + (1 - 2x p (l - x p )) 



In 



1 - z c (t(x p ) )^,_, 
z c (t(x p )) 



— c<t(x p )<l~2c 



(26) 



where is the characteristic function of the set specified in A restricting the integration to 
that set and t(x p ) is the variable t introduced above given by 



t(x p ) 



x B 1 



1 - x B 



(27) 



The logarithms depending on the factorization scale that cancel part of the scale dependence 
of the parton densities are indicated explicitly. There are no such terms depending on the 
renormalisation scale because the Born term is of O{a s ). 

The corresponding terms for a virtual photon with longitudinal polarization arise at 0(a s ). 
They are given by 



d4 1+1)M 

dx B dy 



2N f 

E 



dx r 



=i Jx b 




1 1 4(l-y) + l + (l-y) 2 



2 S B x B 



2y 2 



and 



4 1+lUg 
dx B dy 



2N f l / 

i=l Jx B X P \ X P 



x p - [xp (1 - 2z c (t(x p ))) S c < t(:rp) < 1 _ 2c 



(28) 



2] 2 1 1 4(l-y) + l + (l-y) s 
f l a 2S H x B 2y2 



27T-4 



1 «s 

2 2^ 



2x p (l 



2x p (l - Xp) (1 - 2z c (t(x p ))) S c <t( :rp )<i-2 C 



(29) 



for incoming quarks and gluons, respectively. 



Now the jet recombination scheme ambiguity which arises from the problem of the mapping 
of a "jet phase space" of massive jets in NLO onto a phase space of massless partons in leading 
order is discussed. The real corrections for the (1+1) jet cross section have been obtained by an 
integration of the Born terms for the production of 2 partons over some region of phase space 
specified by the jet cut definition. To cancel the infrared and collinear singularities, one has to 
define effective (1+1) jet variables that allow the identification of the corresponding singularities 
in the virtual corrections and in the contribution from the redefinition of the parton densities. 
This procedure is a map of the (1+1) jet region of the 2-parton phase space onto the (1+1) 
jet phase space. In principle this map is ambiguous, but fortunately the the difference between 
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different maps is for vanishing jet cut. Here the following scheme is tacitly assumed: If a cluster 
algorithm is applied to a 2-parton event that looks like a (1+1) jet event, then the final result of 
this clustering will be two clusters with momenta pa and pb, one of them being massless and one 
of them being massive. The (1+1) jet Born term and the virtual correction always result in two 
massless clusters pc and pp. The identification is done by mapping (pa +Pb) 2 onto (pc +Pd) 2 - 
In this special case both expressions are equal to W 2 . This quantity can be determined from 
electron variables alone, so what is essentially done is the identification of processes in which the 
electrons have the same momenta. This, however, may correspond to very different situations 
on the parton level, since e.g. a radiated final state gluon can be collinear to the remnant jet 
resulting in a remnant with a large momentum or can be collinear to the outgoing quark. The 
situation is similar in the case of the corrections to the (2+1) jet cross section. In that case, 
however, the identification involves two variables and not only one, and so the identification 
using W 2 is not sufficient. 

4 (2+1) Jets: Virtual Corrections 

In this section the calculation of the virtual corrections to the production of (2+1) jets is 
described. The one-loop corrections to the graphs from fig. |5| are given in fig. |(| The one-loop 
diagrams from fig. [?] contribute to the wave function renormalisation and are taken into account 
in the counter term. To obtain the 0(ot 2 ) corrections the diagrams in fig. |6| must be multiplied 
by the Born diagrams. The resulting topologies can be divided into three classes: 

(I) QED-like graphs with colour factor NcCp, 



resulting in terms proportional to the colour factors NcCp and N C C F . The sum of the virtual 
0(a 2 ) corrections averaged over the colour degree of freedom of the incoming parton is 



(II) QED-like graphs with colour factor NcC F \C f — \Nc 



(III) non-abelian graphs with colour factor 



\N 2 C F 



v, q mc. 



L 1 L 2 Q J 2 8(1 - e){C 2 F E 1>q - \n c C f E 2 ^ 




(30) 



v, g mc. 



L X L 2 Q|8(1 - e){ l -C F E^ g - \n c E 2 ^ 




(31) 
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where 



r(l-e) 



-q 2 -ir]l r(l-2e)' 



(32) 



Here Nf is the number of active quark flavours in the fermion loops and fx 2 is the renormalisation 
scale (it is understood that the running a s is always evaluated at the scale fi 2 ). The contributions 
-Ei,Q) E2^q, Ei^g and ^2,g come from the calculation of the matrix elements of the graphs in 
fig. g. The explicit expressions for the £j q and E\ g are collected in Appendix [b[ The trace 
calculations of the matrix elements were done with the help of REDUCE 24]. Then the loop 
integrals were performed by an insertion of one-loop tensor structure integrals (see appendix 
[A]). Here one has to be careful with respect to the imaginary parts of Spence functions and 
logarithms which are important because q 2 < 0. The results obtained here have been checked 



against those of the e + e -case [25]. This was possible because all infinitesimal imaginary parts 
from the propagators were kept in the formulae. 



In (30), ([?!]) the counter terms (in the MS-scheme) to cancel UV singularities (see ||25|| ) 
is already added. Some 1/e 2 - and 1/e-poles remain. These divergences are due to the IR 
singularities of the loop corrections. 



For the processes with an incoming quark the following variables are defined: 

PiPo P2P0 



'-q ■- 



poq 



poq 



a 



(33) 



The divergent parts are 



v, q inc. 



L 1 L 2 Q J 2 8 (i _ e ) Cp T q 



C F 















1 1 




- + - 


fln^- 


e* e 


V x B 




+ 2 




(34) 



For the processes with an incoming gluon the variables are similar: 



PiPo 
poq ' 



P2PQ 

poq ' 



Xg 

a 



(35) 



The divergent parts are 



tvH, 



v, e inc. 



L 1 L 2 Q J 2 8(l-e)-T 9 
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They will cancel against divergent terms from the real corrections. 

5 (2+1) Jets: Final State Real Corrections 

To 0{a 2 s ) one has to consider the contributions from the Born terms in fig. |8] integrated over 
the (2+1) jet phase space region in the 3-parton phase space. 

Again there are graphs with an incoming quark and an incoming gluon. The generic diagrams 
are shown in fig. ^. There are, of course, additional contributions with incoming antiquarks; 
their structure is identical to the quark- initiated processes. 

The integrations become singular if the integrand contains a propagator whose denominator 
vanishes in the integration region. The method of partial fractions is used to separate initial 
and final state singularities. This allows the identification of the terms proportional to O(c ), 
O(lnc) and C(ln 2 c) (c is the jet cut). In this section the final state singularities are considered, 
the initial state singularities are treated in the next section. 

For the process 

e~(k) + proton(P) — > e~(k') + target remnant(p r ) + parton(j>i) + parton(p2) + parton(p3) (37) 

a parametrization of the phase space of the outgoing partons with momenta pi is needed. The 
target remnant with momentum p r = (1 — £)P is described by the variable £. The parametriza- 
tion is chosen such that the integration over the region s±2 < c (px and P2 being collinear or 
pi being soft or p2 being soft) is simple. In close analogy to calculations in e + e~-annihilation 
it is reasonable to describe the two particle phase space of p\ and p2 in the CM frame of these 
momenta (see fig. [l0|). 



Let po be the momentum of the incoming parton. One can define a variable z by 

P0P3 
z := 

poq 



(38) 



that describes the phase space of p% (the azimuthal dependence is contained in the lepton phase 
space, and the remaining third integration is trivial because of energy conservation). One can 
define a polar angle 9 given by 6 := Z(pi,pb) i n the CM frame of p\ and p2 and an azimuthal 
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angle ip by the angle between the planes spanned by po, p\ and po, ps, respectively. Let \ be 
defined by \ '■= ^(PcPs) an d 

b := -(l-cos0) (39) 

d : = i(l-cosx) (40) 



b + d-2bd- 2Jb{\ - b)d(l - d) cos <p. (41) 



With the normalisation factors 

vr4T(l-2e) [* . _ 2e 



r 2 (i- e 

T 2 (l-e) 



N v := — T^TTi = / sin V^P ( 42 ) 



o 



* : = fF^-Jt W1 * (43) 



the 3-parton phase space in d = 4 — 2e dimensions is 

/ dPS(3) = / 128^2- 2e) Sr2EdSl2Z " £ " Xs)(1 " Z) " Sl2) 



128vr 3 r(2 - 2e) 
■dz—— sin _2e tpdip— (6(1 — 6))~ e db 



dPS* (2) /" dPS* (r) . (44) 
dPS*( 2 ) is defined in eq. (O), and 



dPS*W = a^-a)L 2 — m" 2 ^^ a*-^— (45) 

a s 2 87r z rr p ^ 1 — 2e 

= fl~ g 77 S12 V1 rf 2 e dri2-^sin- 2e ^ 

•i-(6(l-6))- e d6, (46) 

(47) 



^pg*(2) gorging the variables z and t = S12/W 2 that are identified with the corresponding 
(2+1) jet variables. The phase space for 3 particles factorizes as a product of a phase space 
for 2 particles dPS*^ r ^ and an effective phase space for a particle and a cluster dPS**- 2 - 1 . The 
invariants can be expressed in the variables t, z, b, x p = x B /£, <p and r±2- 

7*01 = (1 - z)b, 

r 02 = (l-z)(l-6), 

ro3 = z, 

ri3 = (l-Xp- r 12 )e, 

r 23 = (l-x p -r 12 )(l-e). (48) 
14 



The variable d in eq. (|40|) is given by 

d 

The phase space boundaries are 



r\2 



1 — z 1 — x p — ri2 



(49) 



x B E[0,l], S€[xb,1], i€[0,l], 

^g[0,tt], 6 €[0,1], n 2 € [0,(1- x p ) (1-*)]. 



(50) 



The products of the diagrams in fig. || with the complex conjugated diagrams can be classified 
in classes A, B, H with different colour factors (see tab. [I]). An explicit calculation shows that 
the colour classes G and H are regular when integrated over the 2-particle phase space cZPS*^ r ^ 
and therefore vanish for c — ► 0. Therefore these classes are not considered here. 





i 


II 


III 


IV 


V 


VI 


VII 


I 


A 


* 


* 










II 


B 


A 


* 










m 


c 


E 


D 










rv 








F 


* 


* 


* 


V 








G 


F 


* 


* 


VI 








G 


H 


F 


* 


vn 








H 


G 


G 


F 



Class 


Colour Factor 


A 


C 2 F N C 


B 


N C C F (C F - N c /2) 


C 


-1/2N&C F 


D 


N' C C F 


E 


l/2N' c C F 


F 


1/2N C C F 


G 


N C C F (C F - Nc/2) 


H 


l/2N c C F 



Table 1: Colour classes and colour factors of the real corrections. 
The calculation of the spin sum for external gluons has been performed with the formula 



d-l 



(51) 



A=0 



To cancel the contributions from scalar and longitudinal gluons one has to subtract diagrams 
with external ghost lines. The longitudinal and scalar contributions then drop out because of 
the Slavnov- Taylor identities (see (2t| ). 



The matrix elements have been calculated with REDUCE in d = 4 — 
principle they could be obtained from the results in e + e~-scattering [27, 28 



2e dimensions. In 
psfl . However, the 

results of the procedure to obtain partial fractions are different here since some of the invariants 
pick up a sign because of the crossing prescriptions. Here only the results of the factorization 
of the IR divergent terms 



■^singular 



KT, 



q/g- 



(52) 



are stated. T q / g is the Born term 



21| ) with incoming quark and gluon, respectively, and if is a 
singular kernel whose integration is divergent in d = 4 dimensions. In d = 4 — 2e dimensions the 
result of the phase space integral is of the form a/e 2 + B/e+ C. A and B do not depend on the 
invariant mass cut c. C contains terms of the form In c and In 2 c which diverge for c — > 0. The 
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Class 


incoming parton 


product of diagrams 


colour factor 


F x 


quark 


I I, II II, II I 


NoCl/Nc 


F 2 


quark 


h i, m i, m il 


(-1/2)N£C F /N C 


F 3 


quark 


m i, in ii 


{-l/2)N%C F /N c 


Fa 


quark 


iiiiii 


NlC F /Nc 


F 5 


quark 


IV-IV, V-V, VI- VI, VILVII 


{l/2)N c C F /N c 


F 6 


gluon 


LI, II II, II I 


N C C 2 F /(2N C C F ) 


F 7 


gluon 


h i, m i, m il 


(-1/2)N%C F /(2N C C F ) 



Table 2: Colour factors. 



contributions from the final state singularities are divided into seven classes (see tab. ||). The 
results for the traces of the hadronic tensor tiH Fi for the seven colour classes are given explicitly 
in appendix |C|. 

The singular kernels are integrated over the (2+1) jet like region in phase space. The 
condition for partons j, k to form a jet is Sjk < cW 2 . This condition can be rephrased in terms 
of the variables 

'* - ik (53) 

as 

rjk < c — x p . (54) 

The phase space boundary is given by rjk < (1 — x p )(l — z). So the (2+1) jet region is specified 
by ' ^ 

fjk < ol := min |(1 — x p ){\ — z), c —x p \, 6 6 [0,1], <£G[0,7t]. (55) 

^ Xb ' 

The phase space integrals are listed in appendix One obtains 

/ (iPS* (r) trif Fi = a5(£ - a)LiL 2 8(l - e)F t . (56) 

7(2+1) jet 

for the corrections from the final state singularities. The explicit expressions are collected in 
appendix |E[ 



6 (2+1) Jets: Initial State Real Corrections 

In this section the contributions from the initial state singularities are described. In the case of 
the final state singularities the momentum fraction £ of the incoming parton is a fixed parameter 
in the matrix elements. In the case of the initial state singularities, however, there is the 
additional problem that £ is an integration variable and that £ is an argument of the parton 
densities /«(£)• Although there is no expression for in a closed form, the problem can be 
solved |Q] by a Taylor expansion around the singular point and the well known "+" -prescriptions 

n 

D + (g) := f 1 duD(u) (g(u) - g(l)) . (57) 
J o 
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A similar parameterization of the phase space as in the case of the final state singularities is 
used, but the effective (2+1) jet variables are defined in a different way. Let po be the momentum 
of the incoming parton and P3 the momentum of an outgoing parton. For the calculation of the 
singularity resulting from S03 — > p\ and P2 are the momenta of the partons to be identified 
with the effective (2+1) jet momenta. Variables z', z and t are defined by 

, POPS POPl . S12 ,- Q , 

z := , z:= , t:=—». (58 

p q p q W 2 

With the conventions from Section [| one obtains 

6=^, (59) 

A factor of unity 

1 = / " O (60) 

is inserted in eq. (041) to display a change of variables 



1 



l—z' X—Xb 



a := 1 " z . L ~ XB (61) 



1-2' 



explicitly. As a result one obtains 

J dPS (3) = j dPS* {2) J dPS* {r \ (62) 

/dPS*W = 8^-i')L 2 —^\^W 2 H{z')dn u (63) 
J a s 2 8n z 



1-e / / 
/\-2+2e ' « \ / - 



H(z') = (1 - ^"^(l - j— ^-T^J =1 + ^'), (64) 

r(i-&) /i-a*\-« m--,,.^ i „ ; „- 2e 



^ = _^ £ a ^ (1 - ty- e a- £ daz'~ e dz' — sin~ 2e wd<p. (65) 

r 2 (l - e) V 2b J N v 

The ^-function <5(£ — £') can be used to perform the integration over the parton densities, d/sj 
is the measure for the singular integrations. The invariants tjj := Sij/W 2 can be expressed in 
terms of the phase space variables: 

to2 = (y - C)(l - z - z'), 

t03 = (V- C)z, 

t\2 = t, 

hs = (l-C-i)e, 

t 23 = (l-C-t)(l-e), (66) 



where 



1 - x B 



C = (i-^)(i- T 47)' ( 67 ) 



b + d+2bd- 2Jb{l - b)d{l - d) cos ip, (68) 



1 - z' 1 - ^ 



(69) 
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It should be noted that 

£ = a+ (l-a)a + 0(z'). (70) 
The phase space boundaries are given by 

afe€[0,l], ^€[^,1], tE [0,1], ze[0,l], 
<p€[0,ir], a €[0,1], < z < min{l - z,l - t}. (71) 

The factorization of the divergent parts is performed in the form 

M singular = K ' T ?Ar ( 72 ) 

up to terms that vanish for c — > after the integration. Therefore one can, for example, set 
H(z') identically to 1. 

The contributions from the initial state singularities are divided into seven classes (see tab. FJ). 
In the case of the initial state singularities the incoming parton of the (3+1) jet graph is not 
necessarily the incoming parton of the (factorized) (2+1) jet process. This is the familiar fact 
that the quark parton densities modify the evolution of the gluon density and vice versa. In 
tab. |3] a list of the incoming partons of both processes is added. The explicit expressions of the 
singular kernels triJ/ ; are collected in appendix ||]. 



Class 


(3+1) 


(2+1) 


product of diagrams 


colour factor 


h 


quark 


quark 


li, nil, ii i 


N C C F /N C 


h 


quark 


quark 


ii i, in i, mil 


(-1/2)N%C F /N C 


h 


quark 


gluon 


IVdV, V-V, VI- VI, VILVII 


(1/2)N C C F /N C 


h 


gluon 


quark 


LI, ILII, II I 


N C C 2 F /(2N C C F ) 


h 


gluon 


quark 


h i, iii-i, m il 


(-l/2)N 2 C F /(2N c C F ) 


h 


gluon 


gluon 


ii i, in i, m il 


(-l/2)N 2 C F /(2N c C F ) 


h 


gluon 


gluon 


mill 


N 2 C F /(2N C C F ) 



Table 3: Colour factors. 



The boundaries of the (2+1) jet phase space region are given by the invariant mass cut 
condition. For the initial state singularities, the variable z' is the crucial variable that determines 
the singularity structure. The matrix elements become singular for z' = 0. Let "r" be the label 
for the target remnant jet. The invariant t r j is given by t r j = ((1 — £)/£) ioj- If Pr and P3 are 
the momenta combined into a jet, then the invariants of the effective (2+1) jet event are £12, 
tir3, ^2r3- The cut conditions read 

t r3 = ^—^z'<c, (73) 
1 - x B 

tu = t>c, (74) 

1 - £ 

hrs = t 13 + t lr + t 3r = (l-(-t)e + - -{z + z')>c, (75) 

1 - x B 

t 2 r3 = t 23 + t 2r + t 3r = (l-{-t)(l-e) + ^—^-(l-z)>c. (76) 

1 — Xg 
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Especially the last two of these conditions are too complicated to be used in an analytical 
calculation since they involve a restriction of the azimuthal angle integration. Therefore the 
contributions from the initial state singularities are integrated up to the phase space boundary 
(so ip £ [0, tt],z' G [0, min{l — z, 1 — t}]) by keeping the effective (2+1) jet variables z and t fixed. 
The (3+1) jet contribution is then subtracted after a numerical integration. Since the integral 
including the parton densities cannot be performed analytically, this is not a serious restriction. 
The shape of the phase space regions in the (z', <r)-plane is given in fig. |ll]. 

The poles and double poles in e characterising IR (z' = a = 0) and collinear (z' = 0) 
singularities can be calculated in the integration over the full (z' , cr)-plane. The phase space 
integrals are given in appendix^. In the formulae given there the upper limit of the ^'-integration 
is /3 = min{l — z,l — t}. In the integrals a is used as an integration variable. The integral 
involving the parton densities is of the form 

V<«y v --\ q V). (77) 







Here D is a generalized function depending on a and the other jet variables. f {£,((?, z', . . .), Q 2 ) 
is expanded in a Taylor series in z' and all terms of order O(z') that do not contribute in the 
approximation used here are neglected. One obtains 

/(£(<r, z', . . .) = /(a + (1 - a)a, Q 2 ) + O(z'). (78) 

a is the momentum fraction of the incoming parton of the factorized Born term. With the 
definition 

u ■= 1 = >,? \ + °^ ( ?9 ) 

£ a + (1 — a)a 



eq. ( |77| ) can be rewritten as 

f ld -^f^,Q 2 )D(a(u))-^. (80) 
jo u u I — a 

Since D is a generalized function, one has to take care for the boundary terms of the variable 
transformation a — > u. 

Finally one obtains for the real corrections from the initial state singularities 

/ dPS* (r) tr# 7 . = f 1 —i5 - -) £i£ 2 8(l - e)iv (81) 

7(2+1) jet u (3+1) jet 1 Ja u \ ^ u) 1 2K >% 17 

The explicit expressions for the I{ are given in appendix |l|. 



7 (2+1) Jets: Finite Cross Sections 

In the preceeding sections the calculation of the Born terms 0(a s ), the virtual corrections 0(a 2 ) 
and the real corrections 0(a 2 ) has been described. In the sum of the virtual and real corrections 
the IR singularities cancel, and the remaining collinear singularities are absorbed into the parton 
densities by the redefinition eq. (|14|). The final cross sections is then free of divergencies. The 



19 



partonic cross sections must be multiplied with the different flavour factors of the 14 classes of 
diagrams and integrated over the momentum fraction of the incoming parton. Let the charge 
of the quark of flavour i be qi = Qie, where i = 1 stands for d-quarks, i = 2 for d, % = 3 
for u-quarks, and so on. Let /j(£,M|) be the parton density of flavour i, f g (£, M?) the gluon 
density and Nf the number of flavours. Then the flavour factors are 



H Fl , Hp 2 , H Fs , Hp 4 , H Il , Hi 2 : ^ Q 2 fi (£, M 2 ), 

i=l 



Nt 



Hf 6 , Hp 7 , H l4 , H Is , H I(i , H Ir : Qli-l/s(£> Mj), 

N f 2N f 

H h : E^-lE/i(^ M /)" ( 82 ) 

i=l j=i 

Here the factors Qj are included which are already present in the terms Hp i and Hj.. 
The Born terms and virtual corrections with incoming quarks are multiplied by 

^QViitMf), (83) 

i=l 

those with an incoming gluon by 

Nf 

Y.Qli-M^M}). (84) 

i=l 



8 Numerical Results 



In this section the results of the numerical evaluation of the jet cross sections are presented. 
The finite (1+1) jet cross section has been calculated in Section ||[ the (2+1) jet cross section 
in Section [?]. 

The matrix elements for (1+1) and (2+1) jet production are implemented in the program 



PROJET 3.3 [ 30 1 which uses the multidimensional adaptive integration routine VEGAS |31, 
32 [ for the numerical integrations. The parton density parametrizations are from the package 
PAKPDF p3| . PROJET 3.3 allows the integration over bins in x B , y, W 2 and Q 2 . Furthermore, 
acceptance cuts on the angles of the outgoing lepton and the outgoing jets in the laboratory 
frame can be applied. 

To be definite, the following parameters are used. In the case of HERA the CM energy is 
-^CM = 295 GeV, whereas the CM energy of the E665 experiment is -E^m = 31 GeV. In the 
latter case the lepton phase space is restricted by 2 GeV < Q < 5 GeV, 20 GeV < W < 40 GeV 
and 0.05 < y < 0.95. Unless otherwise stated, the renormalisation and the factorization scales 
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are set to /i = Mf = Q, and the parton densities are from the MRS, set D- parametrization 



[34], which is presently favoured by structure function measurements at HERA. The number of 



flavours in the final state is set to 5. 



The dependence of the jet cross sections on the jet cut c is shown in fig. 12. The jet definition 
scheme is Sjj < cW 2 (sij < M 2 means that two clusters have to be combined if their invariant mass 

is smaller than M 2 and have to be considered as two separate clusters if their invariant mass is 
larger than M 2 ). The phase space of the outgoing lepton for HERA is assumed to be 0.001 < 
x B < 1, 10 GeV < W < 295 GeV, and 3.16 GeV < Q < 10 GeV (a), 10 GeV < Q < 31.6 GeV 
(b), 31.6 GeV < Q < 100 GeV (c). The graph in fig. [jj (d) is for E665. The (1+1) jet Born 
cross section does not depend on the jet cut, because the condition s r i > cW 2 is always trivially 
satisfied for c < 1. The (1+1) jet cross section in NLO decreases with decreasing cut c, because 
the total cross section to 0{a s ) is independent of c while the (2+1) jet cross section to 0(a s ) 
strongly increases with decreasing cut. The (2+1) jet cross section in NLO is comparable to the 
(2+1) jet cross section on the Born level as long as the jet cut is not too small. At c ~ 0.007 
the NLO (2+1) jet cross section starts to decrease, and will go to — oo for c — > because of 
dominant terms ~ — In 2 c in the NLO corrections. If one thinks of cW 2 as a new scale in the 
cross section, this behaviour (an extremum of the cross section at some value of this scale) can 
be interpreted as a stabilsation with respect to a change in this scale. If the difference between 
the cross section on the Born level and the cross section in NLO is too large this is a sign of 
the breakdown of fixed order perturbation theory. However, if c > 0.01 this does not seem to 
be the case, and therefore the region of values for c for phenomenological studies should start 
here. The contribution from the longitudinal polarisation of the virtual photon is always small 
(of the order of 20%) compared to the transverse cross section. Therefore using the transverse 
contributions in NLO and the longitudinal contributions on the Born level should be accurate 
(the relative magnitude of the corrections to the longitudinal cross sections are expected to 



be of the same order as in the transverse case, see also [14]). The cut dependence of the jet 
cross sections in all the three regions in Q for HERA studied here is similar, up to the absolute 
normalisation. 

The jet rates #(2+1) m ^S- 13 are defined by #(2+1) = cr (2+l)/ (T tot- I n the case of Born 
terms, cr^ot is given by the total cross section to 0(a®), which is equal to the (1+1) jet cross sec- 
tion on the Born level. In the case of the NLO corrections, cr^ot is given by the total cross section 
to 0(a s ), which is equal to the sum of the (1+1) and (2+1) jet cross sections to 0(a s ). This defi- 
nition has the advantage that the denominator is always independent of the jet cut. #(2+1) Born 
is strongly increasing with decreasing jet cut c. #(2+1) NLO ^ s smaner than #(2+1) Born m 
the bin of smaller values of Q (a). In the bin of larger values of Q (c), #(2+1) NLO * s l ar g er 
than #(2+1) Born ^ or c ^ 0.015 and smaller for c < 0.015. 



Now two different jet definition schemes are compared, fig. 14. Here we use the HERA CM 



energy. The lepton phase space is given by 3.16 GeV < Q < 20 GeV, 10 GeV < W < 295 GeV. 
In (a) the scheme is Sjj ^cW 2 , and in (b) the cut condition is Sjj <c 

(V^v^y 1 "^) with 

a = 0.8, (3 = 0.2. The (2+1) jet cross sections on the Born level in the scheme (b) are always 
larger (for the same value of c) because the absolute scale of the cut in (b) is always smaller. 
Furthermore, the NLO starts to deviate considerably from the Born level at c ~ 0.02 in the 
jet definition scheme (b). Therefore, if this scheme is used, larger values of c are advised. If 
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the parameter /3 is too small, then the NLO cross sections are frequently negative (for small 
absolute values of the cut scale, fixed order perturbation theory breaks down, this is similar to 
the case of small c in the cW 2 scheme, see above). 

The results for the scale dependence of the cross sections at HERA energies are shown in 



fig. 15 (a)-(c) (0.001 < x B < 1, 5 GeV < Q < 100 GeV, 10 GeV < W < 295 GeV, c = 0.02). 



In principle, the renormalisation scale p and the factorization scale Mf are arbitrary. These 
scales give rise to logarithms of the form ln(^ 2 /M 2 ) and ln(M 2 / M 2 ) in the cross section, where 
M is some mass scale in the process. The logarithms are potentially large and spoil perturbation 
theory if the renormalisation and factorization scales are not of the same order of magnitude 
as M. In order to study the behaviour of the cross sections for a change of the scale, the 
renormalisation scale (fig. ^ (a), (c)) and the factorization scale (fig. |l5| (b), (c)) are varied in 
the form of pQ, where p is a parameter in the range 0.2 to 5. If only the renormalisation scale 
is varied (a), the (1+1) jet cross section on the Born level is constant, because it is of 0(oP s ). 
The (2+1) jet cross section on the Born level possesses a large scale dependence of ±40% in the 
range of p given above. The NLO correction to the (2+1) jet cross section reduces this scale 
dependence considerably, because there is a term logarithmic in p that cancels a part of the 
scale dependence of the running coupling in the Born term such that the overall dependence on 
p is (formally) of 0(a s ). The (1+1) jet cross section in NLO is scale dependent because of the 
running coupling constant, and there is no mechanism (i.e., no explicit logarithmic term in p) 
that would cancel this dependence (the reason for this is discussed in Section ||). Fortunately, 
the total cross section to 0(a s ) is less scale dependent than the (1+1) jet cross section. The 
dependence on the factorization scale is shown in fig. [L5| (b). The (1+1) jet Born term is strongly 
scale dependent. Because the parton densities are redefined when the NLO contributions are 
calculated, there is a term that makes the dependence on Mf formally of 0(a 2 ) in NLO. A 
similar cancellation takes place for the (2+1) jet cross section. If both the renormalisation and 
the factorization scale are varied (c), the overall picture is that, compared to the Born level, 
the NLO cross sections are less scale dependent. In fig. |l?] (a) the scale dependence of the 
jet rates is shown.B It is evident that the NLO results have a much smaller scale dependence 
than the results on the Born level. The same graphs for HERA for the range of Q given by 



100 GeV < Q < 200 GeV are shown in fig. 15 (d)-(f), and the corresponding graph for the jet 



rates is fig. ftq (b). The scale dependent cross sections and jet rates for the kinematical region 



of the E665 experiment (for a jet cut c = 0.04) are shown in fig. 15 (g)-(i) and fig. |16| (c). 



Finally the dependence on the parton densities at HERA for two different jet definition 



schemes (fig. 17) is discussed, with parameters 5 GeV < Q < 295 GeV, 10 GeV < W < 295 GeV, 



c = 0.02. The parametrizations HMRS set B pq ], MT set Bl [36] and the more recent ones 
MRS sets DO and D- are chosen for comparison. The two sets of curves in fig. 17 are for the 
cross sections x B da/dx B for (1+1) and (2+1) jets in NLO. The jet definition scheme in (a) 
is Sij^cW 2 . For values of x B smaller than 0.01 the different parametrizations clearly predict 

different (1+1) jet cross sections. The (2+1) jet cross section differential in x B is insensitive to 
a variation of the parametrization. The reason is that in the cW 2 scheme all contributions to 
the (2+1) jet cross section come from £ > c, £ being the momentum fraction of the incoming 
parton, as discussed in Section [| Since c = 0.02, there is very small variation in the (2+1) jet 



J The three curves for the NLO corrected terms should intersect at p = 1. The small difference at p — 1 is due 
to small fluctuations from the spline fit to the Monte Carlo results. 
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cross section because there is not much difference in the parametrizations for £ > 0.02. The 

situation is different in the scheme (b) Sij <c (w a Q^\/S H y 1 a ^ with a = 0.7, (3 = 0.3. Here 

the (2+1) jet cross section receives contributions from the parton densities at £ < c as well, 
and therefore the (2+1) jet cross section depends on the parametrization. Using such a scheme 
might therefore be a possibility to measure the gluon density M 2 ) for small £ via (2+1) jet 
cross sections by a a subtraction of the quark initiated contribution from the total (quark and 
gluon initiated) (2+1) jet cross section. 

9 Summary and Conclusions 

In this paper the calculation of (1+1) and (2+1) jet cross sections in deeply inelastic electron 
proton scattering has been described. The jet definition includes the target remnant jet and is 
based on a modified JADE cluster algorithm. The inclusion of the proton remnant in the jet 
definition scheme is a consistent way to define 'exclusive' jet cross section for the production 
of (n+1) jets because of the possibility of collinear emission of partons in the direction of the 
target remnant jet. 

The cross sections are studied for HERA and E665 energies in detail. The jet cut dependence 
suggests that, if cW 2 is used as the mass scale in the invariant jet definition, the jet cut c 
should be larger than 0.01 to avoid large NLO corrections that could invalidate a fixed order 
perturbative expansion. In the cW 2 scheme, the (2+1) jet cross section depends on the parton 
densities fi(£,M 2 ) for £ > c only, even for very small x B . If one wishes to probe the parton 
densities at smaller values of £, a different jet definition scheme has to be used. In the proposed 
region for c, the NLO corrections are at most of the order of 30%, and even smaller for very 
large Q 2 . 

A set of jet definition schemes that could be useful is given by the cut condition 

Sij <c (w a QP S H y l a ^ , where a and [3 are some parameters in the range of [0,1]. By 

comparing the results for different parametrizations of parton densities it is explicitly shown 
that such a jet definition scheme gives a strong dependence of the (2+1) jet cross section on 
the chosen parametrization. However, it must be studied whether such a jet definition is exper- 
imentally feasible and useful for the determination of the gluon density. 

An important point is the scale dependence of the calculated cross section. It is a general 
phenomenon that leading order cross sections that depend on the strong coupling constant 
and scale dependent parton densities are strongly scale dependent. This leads to a theoretical 
uncertainty because, in principle, the renormalisation and factorization scales are arbitrary 
(although they should be chosen to be of the order of some physical scale in the process) and 
the variation of the cross section with respect to changes in the scales can be interpreted as 
being due to (unknown) higher order corrections because the cross section to all orders must be 
independent of the scales. The NLO corrections usually improve the situation because terms 
arise that cancel part of the scale dependence of the leading order. This desirable feature is 
present in the calculation described here, as has been shown explicitly by a variation of the 
scales as multiples of Q 2 . The scale dependence is significantly reduced. It can be concluded 
that the NLO corrections reduce the theoretical uncertainties of the leading order and should 
provide well defined jet cross sections that could be useful in experimental analyses. 
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A Massless 1-Loop Tensor Structure Integrals 



This appendix contains the results for the tensor structure integrals 



4{0,lM,2^}( r '^) : = 



d d k 



d d k 



(k 2 + ir]){{k — r) 2 + ir/) ' 

{1, /c^, k^k^y 



ho(pi,P2,P3) ■= J d d k 



(k 2 + ir])((k — r) 2 + ir]){{k — r — p) 2 + irj) ' 

1 

(k 2 + ir])((k +P2) 2 + irj)({k - p\) 2 + irj)((k - p\ - P3) 2 + ir]) 



(85) 



needed for the evaluation of the virtual corrections for graphs with massless particles. It is 
assumed that p 2 = p\ = p\ = p\ = 0, but not necessarily r 2 = 0. The integrals are regularized 
by dimensional regularisation (d = 4 — 2e) to take care of UV and IR divergences, ir/ is an 
infinitesimal imaginary part, rj > 0. It is convenient to define the functions 



F(r,p) 
G(r,p) 

H(r,p) 



r 2 + 2pr + ir] \ 



q 2 + irj 



r + ir] 
q 2 + irj 



2pr 



2pr + ir) \ 



l-e 



2pr 



q 2 + ir] I 

I r 2 + 2pr + ir/ 
I g 2 + ir] 



r + i?7 
q 2 + ir] 



l-e 



(86) 
(87) 



2-e 



r + ir; 
q 2 + ir/ 



2-e' 



The calculation with Feynman parametrizations gives (compare p5| ) 



ho{r) 
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r 2 



+ -—F(r,p) + (e + e 2 )G(r,p) 
2pr 

T , x ivr 2 -T(l + e)r 2 (l-e) . 2 . . 
ho(Pi,P2,P3) = K^f. 7T\ (~1 ~ ir l) 



(p^r u + Pu r^+0(e), (93) 



e 2 r(l-2e) v * " (q 2 ) 2 y 12 y 13 
■[l-e(l(y 12 )+l(yi 3 )) 

+e 2 (^ 2 (yi 2 ) + ^(yis) + R(yi2, m)) } + 0(e). (94) 

In I40 the momentum g is denned by q = p\ + p 2 + p 3 , and the invariants y^ are given by 
yij := 2piPj/q 2 . The function R is given by 

R(x, y) = l(x) l(y) - l(x) 1(1 - x) - l(y) 1(1 - y) - S(x) - S(y) + C(2). (95) 

l(x) is the natural logarithm with an additional prescription for arguments on the cut [ — 00, 0] 

l(x) := lim ln(x + sgn(g 2 )sgn(l — x)ir/), (96) 



and S is defined by 
where C 2 is the complex dilogarithm 



S := lim C 2 (x + sgn(q )sgn(l — x)ir)), (97) 



C 2 (z) = - Tdu^-^. (98) 
JO u 

It can easily be seen where the i?7-prescription is important. For q 2 > it fixes the imaginary 
part of the factor (— q 2 — ir])~ e . Expanded up to 0(e 2 ) this gives the well known 7r 2 -terms (in 
combination with l/e 2 -poles) in e + e~-annihilation and in the Drell-Yan process. In deeply 
inelastic scattering (— q 2 — ir])~ e has no imaginary part. However, the functions F, G and H 
give rise to 7r 2 -terms in combination with poles 1/e, 1/e 2 . 



B Virtual Corrections 

In this appendix the explicit expressions for the virtual corrections are given. For the processes 
with an incoming quark the following invariants are defined: 



Zn 1 Zq 

Vqi 

Dig = = ~i 

Xp Xp 

y q9 = 09) 
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where the variables z q and z g are defined in eq. (|33|) and x p = x B /£. One then obtains 
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For the processes with an incoming gluon the following variables are defined: 
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1 - z„ 



(102) 



The variables z q and z q are defined in eq. (|3q). For these processes one obtains 
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The function R' is given by 



R'(x, y) = In |x| In \y\ — In \x\ In |1 — x\ — In |y| In |1 — y\ 
- lim Re (£ 2 (x + ii]) + £ 2 (y + ir/)) + ((2)- 

r)\0 
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C Factorized Integration Kernels for Final State Singularities 



In this appendix the results for the singular kernels from the final state singularities are sum- 
marised. It is convenient to add additional indices to the integration variables r, z and b. Let 



"3 ■- 



PoPj 



where po is the momentum of the incoming parton, and 

Sjk 



r jk := 



(106) 
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If the integration variable b is fixed in the (pj,Pk) CM system, let 



(108) 



Let T q / g (x p , Zj) be the Born term with incoming quark/gluon, expressed in the variables x p = 
x B /£ and Zj. For the singular matrix elements one obtains (including the average over colour 
degrees of freedom for incoming partons and the symmetry factor for identical particles in the 
final state) 
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In this factorization terms that vanish (after the integration) for c — * are neglected. Therefore 
the cut should not be too large in the numerical evaluation. For small values of the cut the cross 
section is dominated by terms ~ ^— In 2 cj . 



D Phase Space Integrals for Final State Singularities 

In this appendix the results for the real corrections of the terms involving final state singularities 
are collected. A measure dfip is defined by 



i, F := ja irr- [l - - j jl _»-(! - *)- £ ^ In- (116) 

Integrals of terms with singularities for r — ► 0, b — * are needed in the integrations over the 
singular region of phase space. iV& and iV"^ are normalisation constants from the phase space, 
and h is an arbitrary parameter which does not show up in the results up to order O(e ). The 
upper limit a of the r-integration has the meaning of a jet cut. The integrals that were needed 
were solved by use of [37, 38, 3£, 40] and are of the following type: 
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The function fi(y) can be checked against the result in [27]. 

E Real Corrections, Final State Singularities 

By means of the basic integrals from appendix [D| one obtains the explicit expressions for the 
final state real corrections: 
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\ l --l + \^~\ + Oiel (127) 
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Fa = 



-CpQ^Tg{x p , Zq) 



1 1/3 



1 — x 



7 3 , a 



+ 7 ( 2~ ln_ ^ ) ~ o ln ~ -lnxpln(l-xp) 



F 7 = 



+ 1 in 2 x p + 1 ln 2 (l - x p ) - In 2 - 2C(2)| + 0(e) 

+ (q < > q), 

- + In x D In ^ + In 2 In 2 



In 

e l-Xr 



1 - Xr 



1 - Zn 



1 Xr 



+i (ln 2 (l - x p ) - ln 2 (l - z q )) + 25 ("j^) } + 0(e) 



+ (q < > q). 



(128) 



(129) 



F Factorized Integration Kernels for Initial State Singularities 

This section contains the singular kernels from the factorization of the matrix elements with ini- 
tial state singularities. Let T q / g (tjk, z%) be the Born term with incoming quark/gluon, expressed 
in the variables tj^ and z\. The singular matrix elements then read 



trH h 



tiH h 



-{l-V-t qg ) 



+2- 



-(l-v-t qg ) + {l-t qg )a / 

Tq(tqg, Zq), 



(1-e) 



trH h = 



-{l-V-t qg )z' + (l-t qg )z q a 

L 1 8 7 r 2 g/x 2e (-16)(l - e)N c C F ^Q) 
1 

-(1 - V - t qg )z' + (1 - t qg )z q <J 
1 - Z a 



(l-V-t qg )z' + (l-tqg)(l-Z q )<J 

a, 

2~7T 



Tq(t qg , Z q ), 



Li87r 2 ^V e Q 2 



^-<^F^Tg(tqq,Z q )(l-e) 



+ 5 1 - < * ) ^ (l + e) 



-(1 - V - tgg) + (1 - t qq )a (l-U-tqq 



-16C F 



1 1 2(1 - v)t qq a(l - t. 



trflV, = 



W 2 Z ' {l- V -tqqYZq{l-Z q 



A _ 2(1 - e) cos 2 1 

Zq) J 



(130) 



(131) 



(132) 
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tvH h 
tvH l6 



tvH 



1 



[1-v-t 



-(1-6) -2- 



99/ 



^(i - tgg) 

(1 - V - t qg ) + (1 - t qg )a 



Tq {tqg j Zq 



+(q <-> g), 
0, 

1 - z„ 



(133) 
(134) 



(1 - 1/ - t 9 g)z' + (1 - %)(1 - Zq)G -(1 - 1/ - tgg)z' + (1 - tqq)z q <J 

'-"-^ Wt q - q ,z q ), (1.35) 



(-(l-l/-ttf) + (l-tfl)(7)' 

16(1 -e)N c ^Tg(t qq , Zq )± 



+ 



(1 - tgg)<7 



-(l-I/-%) + (l-%)(7 (1-U-t 



VI) 



1RAT 1 1 2^(1-^)(1-C-^) / 



e J cos 



(136) 



In these formulae the terms proportional to 1 — 2(1 — e) cos 2 ip are stated explicitly, since for 
invariant mass cuts the integration over the azimuthal angle is (in general) not over the range 
[0, 7r], and as a consequence 

J dip sm~ 2e cp(l - 2(1 - e) cos 2 tp) (137) 

does not vanish.! Additional factors of (1 — e) and 1/(1 — e) in tnHi 3 and triJ/ 4 , respectively, 
are due to the fact that in d dimensions a gluon has (d — 2) (physical) helicity states, whereas a 
quark has only 2, and the incoming particle in the (3+1) jet Born term is not the same as that 
in the factorized (2+1) jet Born term. 



G Phase Space Integrals for Initial State Singularities 

Here the results for the integrals needed for the real corrections with initial state singularities 
are stated. A measure d\ii can be defined by 

/ *« = r ^- r i: % * r ^0) a (^r) ~' (i - t)i_ '- <i38) 

For an arbitrary C°°-function g : [0, 1] — ► C let 

F(a)[g] := [' daF(a)g(a). (139) 
Jo 

6 These terms therefore do not factorize in the familar form of a singular kernel multiplied by the Born term. 
However, this does not affect the factorization of the singularities because the integration region for <p is [0, n] in 
the limit z' — > 0. Therefore the resulting contribution is not divergent although the integrand is singular in this 
limit. 
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Then the following integrals are given by 

i y 



h(y) 



i 4 



h(y) 



k(y) 



a(l-t) 



+ yo- 



2e 2 e V 2 x B J 4 x B 



H-(-f)}« 



+ { --+ln 

e x B 

'ln(l + f)' 



b] 



+ 0(e) 



i 2 := J dfi ig (a)^ 



a(l-t) 



a(l-i) 



^] + (ln(^^)W 
1 1/ 



+ 0(e), 



e V 

J dmg{&)—a 
a(l-i) 



:= /" dmg{o) — 



i y 2 a 



+ 0(e), 



a(l - t) 



1 / y 2 a 



dpigfa) — 



\(l + yaf 

I y 2 



[.9] 



+ 0(e) 



z' (1 + yaf 



= a(l - t) 



e Wl + ya) 5 



V 2 , (l-a)/3a 



In 



(1 + yu) 2 x B 



[5] 



+ 0(e) 



Here the distributions 



<5obl : = 9(0), l[g] ■= C dag(a) 

Jo 



have been used. 

The "+" -prescription for the a- integration is defined by 



f 1 daD + (a)g(a) := C daD(a) (g(a) - g(0)) . 
Jo Jo 
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H Real Corrections, Initial State Singularities 



The explicit expressions of the sum of the (2+1) and (3+1) jet contributions from the initial 
state singularities read 



h = Q 2 3 T q (t qg ,z q ) 



M l 



h 

h + h 



1 o 2 \ 



+c 2 F 



\S{\ - u) + - (- In ^5(1 -u) + h(l - u)) 

e e \ x B 1 J 



h = 



\x B x B (5 
Q]T,(t m ,z,)-(-±N c C F ) 



Ilnl 



QjTg {t qq , Z q ) 



z *xn ( az i Q ( x ~ z <^ z i l ~ Z( > 



Q 



M) 



1 



a(3 

Xg 



+ 0(e), 



QjTq(tqg, Zq) 

f or\' 

+ c F s, ( a ± 

\x B 

0(e), 

QjTg(tqq, Z q ) 



i , Q 2 



+ (q^q) + 0(e), 



q^ g {u)C F 



1 



Q 2 



1 



+ 



r«5(l - u) 



a 2 Z q (l - Z q ) 



+ 1 [ i n ~~«v- o V S (1 - u) - t=- ( ^N C - ^N F ) 5(1 - u) 



x B 



11 



^ ' aZ q q(l - Zg) Zq \-^_z q (l(3 
x B ' x B '/?'/? 



N c V 6 



(148) 



+ 0(e), (149) 



(150) 



(151) 
(152) 



(153) 



The functions Pb^a are the Altarelli-Parisi kernels and the functions Si are defined by 

SM,B,C) = R 2 + R 5 (A)-R 6 (A)+(^ln 2 B + 6((2)")R 1 
+2R 3 (B) - 2R 4 {C), 



(154) 
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S 2 (A,B,C,D) = R 1 (^\n 2 A-Un 2 B^+2(R 3 (A)-R 3 {B)) 

-2{R 4 (C)-R 4 {D)), 
S 3 (A) = Re(A) + 2R 7 (A)-2R 10 , 

S 4 (A) = ±R + ±R 5 (A)-R 8 (A), 

S 6 (A,B,C,D,E) = -R l ^ln 2 A + ^ln 2 B + l2((2) ) j-2(R 3 (A)+R ;i (B)) 

+2 (R 4 (C) + R 4 {D)) + 4R 9 (E) - 4R 6 (E) - 4R 7 (E), 

and the functions Ri are distributions in the variable u given by 



i?0 


= 1, 




Ri 


= 5(1- 


- u), 


R2 


= ] 


- u, 






( In 


K 












1 - 



Ra(X) 
i? 6 (A) 

R 8 (X) 
R 9 (X) 
R10W 



In A 



1 - u 



u 



u 



l-u 
In ( X 



In 1 + A 



l-u 



u 



l-u 



= win A 



u 

l-u 



u 



l-u 



u 



In A 



l-u 



u(l — u) In A 



u 



h In A 



1 - u 



u 



l-u 



u 



(155) 
(156) 

(157) 
(158) 

(159) 
(160) 
(161) 

(162) 

(163) 
(164) 
(165) 
(166) 
(167) 
(168) 
(169) 
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Figure Captions: 



Fig. 1 Diagram for initial state radiation. 
Fig. 2 Diagram for 1-parton production. 

Fig. 3 Diagram for the virtual correction to 1-parton production 

Fig. 4 Generic diagrams for 2-parton production. 

Fig. 5 Born terms for 2-parton production (hadron tensor). 

Fig. 6 Virtual corrections to 2-parton production (hadron tensor). 

Fig. 7 Diagrams contributing to the wave function renormalisation. 

Fig. 8 Born terms of 0{a 2 ) (the roman numbering labels the 7 different colour classes). 

Fig. 9 Generic diagrams for 3-parton production. 

Fig. 10 CM frame of pi and p2- 

Fig. 11 (2+1) and (3+1) jet regions in phase space. 

Fig. 12 Cut dependence of jet cross sections, "ir." stands for transverse, "long." for longitu- 
dinal polarisation of the exchanged virtual photon. Contributions marked as ll NLO" are 
given by the sum of the Born term and the next-to-leading order contribution, (a)-(c) 
HERA kinematics with (a) 3.16 GeV < Q < 10 GeV, (b) 10 GeV < Q < 31.6 GeV, (c) 
31.6 GeV < Q < 100 GeV. (d) is the graph for the E665 kinematics. 

Fig. 13 The ratio cr 2+ i/cr tot . (a)-(c) HERA kinematics with (a) 3.16 GeV < Q < 10 GeV, (b) 
10 GeV < Q < 31.6 GeV, (c) 31.6 GeV < Q < 100 GeV, (d) E665 kinematics. 

Fig. 14 Dependence of the cross section on the jet definition scheme. The jet definition scheme 
is (a) Sij^cW 2 , (b) Si j ^c{w a Q |3 ^/S H ~y 1 ' a '^ 2 with a = 0.8, /3 = 0.2. 

Fig. 15 Dependence on the renormalisation and factorization scales, (a) renormalisation scale 
p = pQ, factorization scale Mf = Q; (b) p = Q, Mf = pQ; (c) p = pQ, Mf = pQ. (d)-(f) 
are the same, but for much larger Q. (g)-(i) are the graphs for E665 energies. To avoid 
scales that are too low for perturbation theory to be valid, the scales are clipped at a lower 
bound of 2 GeV (HERA) and 1 GeV (E665). 

Fig. 16 Dependence of the jet rate o'2+l/ c7 tot on the renormalisation and factorization scales. 
The indices (a)-(i) refer to the indices of Fig. 15 (a)-(i), respectively. 

Fig. 17 Dependence on the parton densities. The plotted cross section is x B da/dx B , the jet 
definition scheme is (a) s^cW 2 , (b) (w a Q? ^/S H ~y l ~ a ~^ 2 with a = 0.7, = 0.3. 
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Figure 3: 
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Figure 4: 
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Figure 6: 
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Figure 7: 
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